full.binomial <- function(n,pi){
a <- matrix(0:n,n+1,1)
b <- dbinom(a,n,pi) 
c <- pbinom(a,n,pi)
result<-cbind(a,b,c,1-c)
rownames(result) <- rep("",n+1)
colnames(result)<-c('y','Pr(Y = y)','Pr(Y <= y)','Pr(Y > y)')
return(result)
}